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ABSTRACT 

With the present paper we conclude the presentation of a semianalytical model of hierarchical 
clustering of bound virialized objects formed by gravitational instability from a random Gaussian 
field of density fluctuations. In paper I, we introduced the basic tool, the so-called "confluent 
system" formalism, able to follow the evolution of bound virialized objects in the peak model. 
This was applied to derive the mass function of objects. In the present paper, we calculate other 
important quantities characterizing the growth history of objects. This model is compared with a 
similar one obtained by Lacey & Cole (1993) following the Press & Schechter (1974) approach. The 
interest of the new modeling presented here is twofold: 1) it is formally better justified as far as 
peaks are more reasonable seeds of objects than the indetermined regions in the Press & Schechter 
prescription, and 2) it distinguishes between merger and accretion enabling us to unambiguously 
define the formation and destruction of objects and to estimate growth rates and characteristic 
times not available in the previous approach. 

Subject headings: cosmology: theory - galaxies: clustering - galaxies: formation 



1. INTRODUCTION 

Observational data on the distant universe are in- 
creasingly numerous and detailed. Comparison with 
the nearby universe harbors important information on 
the formation and evolution of cosmic objects. How- 
ever, our capability of extracting accurate cosmolog- 
ical implications from these observations is severely 
limited by the lack of an analytical model of gravi- 
tational clustering. This is a fundamental lack only 
partially replaced by numerical simulations. 

In the most studied scenario of structure formation 
via gravitational instability from a primordial ran- 
dom Gaussian field of density fluctuations, Press & 
Schechter (1974; PS) derived a practical analytical es- 
timate for the number density of bound virialized ob- 
jects (or more exactly, steady relaxed haloes) of mass 
M at any given epoch. This mass function has been 
shown to agree with A-body simulations (Nolthcnius 
& White 1987; Efstathiou et al. 1988; Efstathiou 
k Rees 1988; Carlberg & Couchman 1989; White et 
al. 1993; Bahcall & Cen 1993; Lacey & Cole 1994). 
But the mass function of objects does not provide all 
the information required in many cosmological prob- 
lems. The rates at which objects grow and the char- 
acteristic times of this process are also needed (e.g., 
Toth & Ostriker 1992; Richstone, Loeb, & Turner 
1992; Kauffman, White & Guiderdoni 1993; Lacey 



& Cole 1993; Gonzalez-Casado, Mamon, & Salvador- 
Sole 1994; Kauffman 1994; Gonzalez-Casado et al. 
1996). Richstone, Loeb, & Turner (1992) used the 
time evolution of the PS mass function to estimate 
the formation rate of objects of mass M at different 
epochs. However, this is not a very accurate estimate 
since the time derivative of the mass function is equal 
to the rate at which objects reach mass M minus the 
rate at which they leave this state, both terms having 
comparable values. 

Following the PS original prescription or a better 
sound version of it using the excursion set formalism, 
Bower (1991) and Bond et al. (1991) derived the mass 
function of objects of a given mass at a given epoch 
subject to the condition of being part of another ob- 
ject with a larger mass at a later time. This condi- 
tional mass function was used by Lacey & Cole (1993, 
LCa; see also Kauffmann & White 1993) to infer self- 
consistent estimates of the instantaneous merger rate 
and the typical age and survival time of objects. This 
clustering model has been shown by Lacey & Cole 
1994 (LCb) to agree with A-body simulations. How- 
ever, as recognized by these authors, there is the for- 
mal caveat that the PS approach on which it is based 
is heuristic. In particular, the seeds of objects are in- 
determined, possibly even unconnected, regions. As 
a consequence, one cannot really count objects with 
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a given mass, but just calculate the probability that 
a given point is in an object of that mass, which is 
at the base of a slight inconsistency in the analytical 
estimate of the formation time for power spectra with 
index n > (see LCa). 

But there is perhaps a more serious problem with 
the PS approach. This deals with one unique process 
of instantaneous mass increase of objects, their cap- 
ture of other more or less massive ones, process which 
is generically called merger. Although strictly correct, 
this viewpoint is at the base of an important short- 
coming. For usual power spectra of the initial density 
field, the number density of objects diverges at small 
masses. Consequently, objects continuously experi- 
ence captures, there being no specific event in this ap- 
proach marking the beginning or the end of any stable 
entity. The age and survival time of an object must 
then be defined rather arbitrarily in terms of some 
relative mass variation along the series of embedded 
objects connecting with it. Moreover, these defini- 
tions are ambiguous since they do not account for the 
fact that the same mass increase can be achieved in a 
smooth quasi-continuous manner or essentially in one 
unique capture at some intermediate time. 

This is particularly disturbing since the notions 
themselves of the formation and destruction of ob- 
jects have been introduced with the aim to specify 
when the instantaneous mass increase of a system 
is important enough to appreciably affect its struc- 
tural properties. In this sense, the distinction be- 
tween "tiny" and "notable" captures would be very 
convenient. Kitayama & Suto (1996) have recently 
proposed a modified version of the LCa model in this 
spirit. However, the PS approach is not well adapted 
to such a distinction which makes the model of Ki- 
tayama & Suto be somewhat inconsistent (see § 2 for 
details) . But this does not imply, of course, that such 
a distinction cannot be self-consistcntly achieved in 
any other approach. 

In tiny captures, the internal kinetic and potential 
energies of the capturing system are not appreciably 
modified from their equilibrium values. So the system 
remains essentially unaltered despite the slight extra 
mass added. In notable captures (possibly the same 
preceding events but considered from the viewpoint of 
the partner object), there is on the contrary a marked 
departure in the internal kinetic and/or potential en- 
ergies of the capturing system. This causes an im- 
portant rearrangement of the system as it relaxes ap- 
proaching the virial relation. We can therefore natu- 



rally say that, in the former case, the final and initial 
relaxed systems are the same object (or equivalently, 
that the object survives to the capture) while, in the 
latter, they constitute two distinct objects. In other 
words, we can define accretion as those tiny captures 
contributing to the (quasi)continuous mass increase of 
objects, and mergers as those notable ones destroying 
and giving them rise. Objects are then well-defined as 
those stable entities surviving, while smoothly accret- 
ing mass, between two mergers marking their forma- 
tion and destruction, and the age and survival time of 
objects would have an obvious unambiguous meaning. 

The problem with this viewpoint is, of course, 
where to place the division between tiny and notable 
captures. Establishing whether there is a significant 
rearrangement of the system as a consequence of a 
capture is equivalent to establishing whether the cap- 
ture causes to it a significant departure from steadi- 
ness. So the difficulty in making that division is 
similar and, in fact, tightly connected with that en- 
countered in making the distinction between relaxed 
and non-relaxed objects (or in identifying the exact 
extent of the former ones): deciding whether a re- 
gion is steady or not. As far as objects are continu- 
ously capturing other more or less massive ones strict 
steadiness is never achieved. Consequently, making in 
practice any of these decisions requires choosing some 
confidence level for assessing the significance of the 
departure from steadiness, which introduces some ar- 
bitrariness. However, the advantage of working with 
theoretical models is that they do not deal with the 
full complexity of the real universe but with some 
suitable idealization of it. For example, the model- 
ing of gravitational clustering presumes it is possible 
to tell, in some unspecified manner, between relaxed 
and non-relaxed systems. So why could it not also 
presume it is possible to tell, in a similar unspecified 
manner, between accretion and merger in the sense 
above? 

As shown in paper I (Manrique & Salvador-Sole 
1995; see also Salvador-Sole & Manrique 1995), mak- 
ing these two sole assumptions in the peak model 
framework leads to one unique, well justified, self- 
consistent, clustering model. The basic tool is a new 
formalism, the so-called confluent system of peak tra- 
jectories, able to follow the filtering evolution of peaks 
tracing the clustering evolution of objects (Salvador- 
Sole & Manrique 1994). This formalism was applied, 
in paper I, to derive the mass function of objects. 
In the present paper, we apply it to calculate their 
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growth rates and times. This clustering model is sim- 
ilar to that developed by LCa in the sense that these 
quantities are derived from the statistics of the contin- 
uous random density field encountered at a fixed ar- 
bitrary epoch after recombination when fluctuations 
are still linear (with the growing factor only a func- 
tion of time at all interesting scales) and Gaussian dis- 
tributed, and the dynamics of dissipationless collapse 
are approximated by the spherical model. However, it 
is physically better motivated since the assumed seeds 
of objects are peaks instead of the indetermined re- 
gions in the PS approach. In addition, it allows us to 
unambiguously define the formation and destruction 
of objects and calculate their characteristic growth 
rates and times. It is important to remark that there 
is no freedom in the model; there are only one or two 
parameters (depending on the power spectrum) gov- 
erning the dynamics of collapse which can be adjusted 
by fitting the real mass function of objects. In par- 
ticular, there is no arbitrariness associated with the 
distinction between merger and accretion. 

But what about the connection between this ide- 
alized model and real (simulated or observational) 
data? As mentioned, the practical distinctions be- 
tween relaxed and non-relaxed and between tiny and 
notable captures require the choice of two, likely con- 
nected, confidence levels. Now, if the model is good 
enough it will agree with empirical data for some spe- 
cific values or some narrow ranges of values of these 
confidence levels or any equivalent pair of parame- 
ters, such as the total mean ovcrdensity of the sys- 
tem, (p)/p with p the mean cosmic density, and the 
relative mass increase of the system AM/M (both in 
principle functions of M and t), for the distinctions 
between relaxed and non-relaxed objects and between 
merger and accretion, respectively. These would be 
the unspecified manners of carrying out those distinc- 
tions implicitly assumed in the model. This is the 
scheme followed by the PS approach for dealing with 
the non-trivial problem of defining and identifying re- 
laxed objects. The model here proposed is intended 
to similarly deal with the problem of defining and 
identifying the formation and destruction of relaxed 
objects. 

The predicted mass function of relaxed objects 
agrees with the PS one recovering the results of N- 
body simulations for relaxed haloes identified as those 
regions with overdensity (p)/p greater than a univer- 
sal threshold (for ft = 1) equal to 178. Cole & Lacey 
(1996) have shown, indeed, that this corresponds to 



the maximum extent of regions in steady state at 
some reasonable (unspecified) confidence level. The 
validity of our model concerning the predicted growth 
rates and times is harder to check because no practical 
distinction has ever been attempted to implement in 
numerical simulations between accretion and merger. 
Nonetheless, a first hint can be provided by the com- 
parison of our results with those obtained from the 
LCa model also checked to agree with numerical sim- 
ulations. 

In § 2 we remind the basic lines of the confluent 
system formalism. In § 3 we derive the instantaneous 
formation, destruction, and mass accretion rates of 
objects. Their typical ages and survival times are 
calculated in § 4. Our results are discussed and sum- 
marized in § 5. The notation used in this paper is 
the same as in paper I. Since many calculations are 
based on results obtained by Bardeen et al. (1986; 
BBKS) we have kept as close as possible to the nota- 
tion introduced by these authors. The main difference 
comes from the fact that we deal with three different 
kinds of densities, while BBKS only dealt with one. 
Especial caution must be made in not mixing them 
up. Firstly, there are the normal and conditional den- 
sity functions of peaks at a fixed filtering scale R per 
infinitesimal ranges of v, i.e., the density contrast S 
scaled to the rms value ao(-R), and other possible vari- 
ables. These were already defined in BBKS and are 
denoted by a caligraphic capital n just as in that pa- 
per. A minor difference with the notation used in 
BBKS is that we specify the fixed value of the filter- 
ing scale R as one parameter. For example, we write 
N(v, R) dv instead of N{v) dv. Second, there are the 
normal and conditional density functions of peaks at a 
fixed density contrast 5 per infinitesimal ranges of the 
filtering scale R and other possible variables. These 
density functions, already introduced in paper I, are 
denoted by a roman capital n with the fixed value of S 
as one parameter. We write, for example, N(R, S) dR. 
Finally, and this is a novelty, there are also density 
functions per infinitesimal ranges of both R and 5 (or 
the corresponding mass M and time t, respectively) 
and any other extra variable. These are denoted by 
a capital n in boldface, for example, N(R,5) dRdS. 
These different symbols, J\f, N, and N, are usually 
accompanyed by one supcrindcx specifying the char- 
acteristic property (if any) of the peaks involved and 
the subindex pk or no subindex at all depending on 
whether all peaks with that characteristic property 
are included or just those of them tracing bound viri- 
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alized objects (see § 2), respectively. Hereafter, we 
assume comoving units. 

2. THE CONFLUENT SYSTEM FORMALISM 

By assumption, a merger is any discontinuous mass 
increase along the temporal series of relaxed systems 
subtending at each step the preceding one. (In prac- 
tice, there must be "some significant discontinuity 
in mass" , or equivalently "some significant structural 
rearrangement" , or still "some significant departure 
from steadiness" between two consecutive relaxed sys- 
tems in the discrete, fine enough, temporal series one 
can consider.) When an object merges we say that 
it is destroyed. Note that if an object is captured by 
any more massive it automatically merges regardless 
of whether or not this more massive object survives to 
the event. If it does not we say that the final system 
subtending the mass of the merging objects is a new 
object that has just formed. 

In contrast, accretion is any continuous (we assume 
it also derivable) mass increase along the temporal 
series of relaxed systems subtending at each step the 
preceding one. (In practice, if there is "no significant 
discontinuity in mass", or "no significant rearrange- 
ment" , or still "no significant departure from steadi- 
ness" along the series of discrete fine steps considered, 
one can regard the mass increase as achieved in a con- 
tinuous and derivable manner.) We say that an ob- 
ject survives as long as it evolves by accretion. When 
an object evolving by accretion is captured by any 
more massive object it merges, hence it is destroyed. 
However, an object evolving by accretion can capture 
less massive objects without being destroyed. In fact, 
its (continuous) mass increase is made at the expend 
of very tiny objects. Only if they capture smaller 
although notable objects (yielding a significant dis- 
continuity in its mass increase) they also merge and 
are destroyed. This is the kind of mergers mentioned 
above giving rise to the formation of new objects. 

Notice that these definitions of accretion and mer- 
ger refer to the object whose evolution is being fol- 
lowed, not to the event itself. This latter can be re- 
garded either as merger or as accretion depending on 
the object under consideration. The inconsistency of 
the model proposed by Kitayama & Suto (1996) arises 
from not accounting for this fact. These authors take 
the model of LCa and make the extra assumption that 
any object which results from another one increasing 
its mass by less than a factor 2 can be "identified" 
to it. This kind of capture is therefore regarded as 



accretion, while that in which the mass increases by 
more than a factor 2 is regarded as a true merger. 
The instantaneous formation/destruction rate of ob- 
jects of mass M is then defined as the intantaneous 
total merger rate for final/initial objects of that mass 
including only those captures of objects less/more 
massive than half/ two-times the mass M. Now, ac- 
cording to the definitions adopted for accretion and 
merger, all these captured objects merge indeed in the 
event, while those more/less massive than half/two- 
times the mass M only accrete. However, not all 
those mergers cause the formation of new objects of 
mass M as presumed; a large fraction necessarily con- 
tributes to the mass increase by accretion of those 
objects identified to the final ones with mass M. 

Coming back to our model, the peak model ansatz 
states that there is a correspondence between peaks 
of fixed linear overdensity in the filtered density field 
at some arbitrary epoch ij and relaxed objects at the 
time t. The overdensity S c is assumed to be a de- 
creasing function S c (t) of the collapse time, and the 
filtering scale R an increasing function R(M) of the 
mass of the resulting objects. Given the assumed dis- 
tinction above between accretion and merger for real 
relaxed objects, that correspondence automatically 
leads to the natural identification and distinction be- 
tween each other of accretion and merger events in the 
filtering process. A peak on scale R + AR, with Ai? 
positive and arbitrarily small, is the result of the evo- 
lution by accretion of a peak on scale R provided only 
that the volume (mass) subtended by the latter is es- 
sentially embedded within the one subtended by the 
former. (Strictly, we should talk about volumes sub- 
tended by "collapsing clouds associated with peaks" , 
i.e., the regions surrounding each peak which enclose a 
total mass equal to that of the corresponding final ob- 
ject at t. But, for simplicity, we will hereafter say the 
volumes subtended by peaks. Likewise, we will say 
"nested peaks" instead of "peaks with nested associ- 
ated collapsing clouds".) Whenever the identification 
between couples of peaks on contiguous scales is not 
possible, that is, whenever there is a discrete jump in 
scale between two consecutive embedded peaks, we 
are in the presence of a merger. It is important to 
mention that the filtering of the density field must 
be carried out with a Gaussian window for the den- 
sity contrast 5 of peaks to diminish with increasing 
filtering scale R as required by consistency with the 
growth in time of the mass of objects. Hereafter we 
adopt this particular window. 
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From this identification among peaks on different 
scales we have that when an object evolves by ac- 
cretion its associated evolving peak traces a contin- 
uous and derivable trajectory 5(R) in the S vs. R 
diagram. (Note that despite the fact that the den- 
sity field is assumed continuous and infinitely deriv- 
able, peak trajectories are only required to be one 
time derivable as the assumed mass increase in time 
of accreting objects.) In contrast, when an object 
merges the evolving peak tracing it becomes nested 
on a larger scale peak with identical <5, which yields 
a discrete horizontal jump in scale of the associated 
peak trajectory in the 5 vs. R diagram. Therefore, to 
compute the density of objects at t in an infinitesimal 
range of masses dM we must calculate the density 
of non-nested peaks with fixed value of 5 appropri- 
ate to t on scales in the infinitesimal range dR cor- 
responding to dM. This density, N(R,S) dR, can be 
obtained from the density of peaks satisfying iden- 
tical constraints although disregarding whether they 
are nested or not, N p k(R,S) dR, and the density of 
these same peaks subject to the condition of being lo- 
cated in a background with the same density contrast 
6 on a different filtering scale R' , N p k(R,5\R' ,5) dR. 
Indeed, these three quantities are related through 

N(R, 5) = N pk {R,6) 

dR'^^-N(R',S)N pk {R,S\R',S), (1) 

IR P 

with M{R) the inverse relation to R{M). Notice that 
this relation is approximated for it assumes the whole 
cloud associated with the peak on scale R' with uni- 
form smoothed density contrast equal to its central 
value. 

Equation (1) is a Volterra type integral equation 
of the second kind for the unknown function N(R, 5). 
This can be obtained by means of the standard it- 
erative process from the known functions N pk (R,5) 
and N pk (R,S\R',5), respectively equal to N pk {v,R) 

(x) 02 R/o-q and N pk {v, R\v' , R') (x) <j 2 R/<to in terms 
of the normal and conditional density functions of 
peaks at a fixed scale R per infinitesimal range of 
v = 5 /ao, calculated by BBKS. In this latter expres- 
sions o~i are the i-th order spectral moments, which 
only depend on R, and (x) and (x) are some averages 
of x = —((72 R)~ 1 dn,o~, which depend on R and 6 as 
well as on R' in the latter case, defined in paper I. 
The mass function of objects at t is then simply 



f 

JR 



with the dependence on M and t on the right hand 
side given by the functions R(M) and 5 c (t). As shown 
in paper I, very general consistency arguments allow 
one to determine the shape of these two functions as 
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(3) 



(4) 



N(M, t) — N(R, S c ) Jg, 



(2) 



with a the cosmic scale factor and q and S c o two con- 
stants governing the exact dynamics of collapse. For 
CDM and the n = — 2 power-law spectra, a good 
fit is obtained to the PS mass function for q equal 
to <~ 1.45 in both cases (in the power law case, the 
larger n, the smaller q) and S c o equal to ~ 6.4 and 
8.4, respectively. (As mentioned in paper I, there 
is no degeneracy in the present formalism between 
these two parameters. It is therefore not surprising 
that a value of q appreciably different from unity is 
coupled to a value of 5 c q so different from the stan- 
dard value of 1.686.) Strictly, equations (3) and (4) 
are only valid for the scale- free case, i.e., an Einstein- 
de Sitter universe (Cl = 1, A = 0) and a power-law 
spectrum. However, as shown in paper I, they are 
good effective relations for other power spectra such 
as the CDM one. Furthermore, following the same 
strategy as in LCa, it should be possible to extend 
the applicability of the model to the cases 0^1 
and/or A ^ 0. We want to emphasize that the val- 
ues of parameters q and <5 c o quoted above, and used 
throughout the present paper to illustrate the practi- 
cal implementation of the model, correspond to those 
giving acceptable fits to the original PS mass func- 
tion (for the standard top hat window and critical 
threshold density) for the same cosmogonies. This 
fitting tends to privilege the small mass end while the 
model will be finally applied rather to massive ob- 
jects. Thus, finer values of these parameters should 
be inferred by directly fitting A^-body simulations in 
the relevant scales. 

Given the characterization of accretion and merger, 
it is clear that the density of objects being destroyed 
(because merging) in a given interval of time is given 
by the density of non-nested peaks which become 
nested along the corresponding decrement in 5. The 
density of forming objects is not so simple to obtain. 
We first need to characterize those filtering events 
which contribute with the appearance of new peaks. 
As just mentioned, when an object merges the non- 
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nested peak tracing its evolution in the 5 vs. R di- 
agram experiences a discrete jump in the scale. But 
this does not mean, of course, that every non-nested 
peak partaking of the same event necessarily experi- 
ences this kind of jump. The largest scale peak will 
the most often just keep on evolving in a continuous 
manner. (This reflects the well-known fact, in gravi- 
tational clustering, that a merger from the viewpoint 
of one given object can be seen as accretion for the 
most massive partner.) But it will sometimes happen 
that the largest scale peak also experiences a discrete 
jump. Then, it will not be possible to identify the 
final non- nested peak with any of its ancestors. This 
appearance of a new non-nested peak therefore traces 
the formation of a new bound object. Hence, the den- 
sity of objects forming (because of some mergers) in 
an interval of time should be given by the density of 
peaks appearing in the sense above along the corre- 
sponding decrement in 8. 

Actually, things are a little more complicated than 
this. The continuous trajectory attached to an ac- 
creting non-nested peak can be suddenly truncated 
(it is not possible to identify the peak at the current 
S with any peak on an infinitesimally larger scale) 
without becoming nested. In addition, peaks can not 
only become nested into larger scale ones but they 
can also sporadically leave their hosts. According to 
the correspondence between relaxed objects and non- 
nested peaks, the preceding events can only be in- 
terpreted as tracing a relaxed object breaking up in 
small pieces and a relaxed object spallating off an- 
other object, respectively. Thus, we are concerned 
with 5 different filtering processes: accretion (contin- 
uous non- nested peak trajectories), aggregation (con- 
tinuous non-nested peak trajectories becoming nested 
into old or newly appeared continuous non-nested 
peak trajectories), fusion (appearance of continuous 
non-nested peak trajectories), and the inverse pro- 
cesses of the latter two, namely, spallation (contin- 
uous nested peak trajectories becoming non-nested 
from surviving or just disappeared continuous non- 
nested peak trajectories), and fission (disappearance 
of continuous non- nested peak trajectories). Since 
spallation and fission are clearly unrealistic processes 
from the gravitational viewpoint, the confluent sys- 
tem formalism cannot be used to follow the clustering 
of individual objects. However, our purpose is not to 
follow the evolution of individual objects, but to own 
a good statistical description of the general clustering 
process. 



Therefore, what exactly traces, statistically, the 
amount of objects that are destroyed (because merg- 
ing) or that form (as the result of some mergers) along 
some interval of time At is the net decrease in old 
non-nested peak trajectories (given by the amount 
of aggregations minus spallations) and the net in- 
crease in new non-nested peak trajectories (given by 
the amount of fusions minus fissions) yielded along 
the corresponding — AS, respectively. Hereafter, we 
will simply refer to these quantities as the net den- 
sity of peaks becoming nested and the net density of 
non-nested appearing peaks, respectively. To calcu- 
late the former we must compute the density of peaks 
at S which in continuously evolving (without fission- 
ing) to S — AS arc found to be nested, minus the 
density of peaks at (5 which are already nested and do 
not fission in the next —AS. Indeed, this latter term 
includes peaks which can end up, at S — AS, being 
nested or non-nested. Hence, that difference really 
gives the density of non-nested peaks at S which have 
just become nested through aggregations, minus the 
density of nested peaks at S which have just become 
non-nested through spallations, as wanted. (We have 
required all peaks at S not to be fissionable because 
aggregation and spallation yield, by definition, con- 
tinuous peak trajectories.) Similarly, to calculate the 
latter net density above we must compute the den- 
sity of non-nested peaks at <5 which do not fission in 
the next — AS, minus the density of peaks at 8 which 
have continuously evolved (without fissioning) from 
peaks at S + AS. Indeed, this latter term includes 
peaks which at 8 become fissionable or not. Hence, 
that difference really gives the density of non-nested 
peaks at <5 that have just appeared through fusions, 
minus the density of non-nested peaks at 8 that are 
just going to disappear (or, in the limit of vanishing 
AS, have just disappeared before reaching 8) through 
fissions, as wanted. (We have required all the peaks 
at 8 to be non-nested because fusion and fission only 
concern, by definition, peaks which are non-nested in 
appearing or disappearing.) Consequently, what we 
only need for the model to be statistically acceptable 
is those two net quantities tracing the densities of ob- 
jects merging and forming along some interval of time 
be positive. 

Let us finally note that mass is usually not con- 
served in individual filtering mergers (more exactly, 
in aggregations or spallations). Indeed, the mass of 
an aggregated (spalled) peak does not yield any in- 
stantaneous finite change in the mass of its partner. 
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Only in the case of fusions (fissions) will the mass of 
the final (initial) peak be close to that of all aggregat- 
ing (spallating) peaks. In any event, the total mass 
of peaks partaking of the filtering processes tracing 
the mergers of objects is certainly conserved statisti- 
cally, as needed. Indeed, the correct normalization at 
every <5 of the scale function of non-nested peaks (see 
paper I) and the fact that the total available mass in 
any representative (large enough) comoving volume 
is fixed guarantee that the the mass shared by peaks 
in that volume is preserved under any instantaneous 
rearrangement through aggregations and spallations 
and the fusions and fissions they experience. 

3. GROWTH RATES 

As just mentioned, the density of objects that 
are destroyed (because merging) in the interval dt is 
traced by the net density of peaks becoming nested in 
the corresponding range — dS. To be more exact, we 
are interested in calculating the net density of peaks 
with 5 C on scales between R and R + dR which be- 
come nested into (i.e., aggregate into, minus spallate 
off) non-nested peaks with 5 C — d5 on scales between 
R' and R' + dR', N d (R -» R' , S c ) dRdR' dS. This 
net density is derived in Appendix B. Now, by divid- 
ing it by N(R, 8 C ) dR we obtain the net conditional 
probability that a non-nested peak with 8 on scale R 
becomes nested into a non-nested peak with 5 C — dS 
on scales between R 1 and R' + dR'. And from this 
conditional probability we can readily infer the in- 
stantaneous destruction (or true merger) rate at t for 
objects of mass M per specific infinitesimal range of 
mass M' (M < M') of the resulting objects, 



r d (M M',t) = 



N d (R^R',5 c ) dR' 



N(R, S c ) dM> 



dS c 
~~dt 



(5) 



with R, R' , and 5 C on the right hand side written in 
terms of M, M' , and t, respectively, through equa- 
tions (3) and (4). (Note that this destruction rate 
involves any possible merger that an object can ex- 
perience with other objects, disregarding whether or 
not this yields the formation of a new object.) 

This destruction (or true merger) rate is plotted in 
Figure 1 for different masses M of the initial object 
and t equal to the present time. As can be seen, it 
shows a very different behavior from the merger rate 
obtained by LCa: the former vanishes while LCa's di- 
verges at small AM/M. This marked difference sim- 
ply reflects the different viewpoints adopted in the 
two models. In the PS approach, any mass increase 



of objects in any given interval of time is necessarily 
due to (generic) mergers, while in ours it can be due to 
mergers or accretion. So both merger rates cannot be 
identical. For small intervals of time, large mass in- 
crements are predominantly caused by mergers rather 
than accretion. Therefore both models should yield 
similar merger rates at very large AM/M as found. 
Conversely, very small mass increments should tend 
to become accretion. In fact, we would even ex- 
pect they all become accretion for AM/M below any 
effective threshold marking the border line between 
tiny and notable captures. So it is well understood 
that our merger rate vanishes at very small AM/M 
while LCa's diverges there because of the divergent 
number of infinitesimal mass captures. (Note that 
what LCb compared with iV-body simulations is not 
their merger rate, but the conditional probability that 
given an object of fixed mass at some initial epoch it 
is incorporated into a larger mass object at a later 
epoch. Since this quantity does not depend on the 
particular definition adopted for merger and accre- 
tion, the good agreement found by these authors be- 
tween theory and simulations does not favor LCa's 
viewpoint as compared to ours.) 

Our merger rate therefore shows the good expected 
behavior at both mass ends. Moreover, the sharp cut- 
off shown at small AM/M points to a division be- 
tween accretion and merger roughly consistent with 
a fixed effective threshold in AM/M (at least for 
given M and t). However, at relatively large values 
of AM/M, our merger rate is significantly lower than 
LCa's which is hard to interpret in terms of our dis- 
tinction between accretion and merger. In fact, this 
effect becomes increasingly marked as we diminish 
M: a central dip develops finally reaching negative 
values. So not only does the correct behavior of the 
model seem to require large M but there is a strict 
lower bound in M for its general physical consistency. 
(Negative merger rates imply more spallations than 
aggregations which is unrealistic in hierachical clus- 
tering.) The value of this lower bound in M is not 
dramatic by itself. At present, it is of the order of 
the mass of small dwarf galaxies (~ 2.5 x 10 9 M 
and - 2.5 x 10 8 M Q , for the CDM and n = -2 
power-law spectra here analized, respectively) and it 
increases with increasing time so that it encompasses 
all relevant cosmological scales (~ 1.8 x 10 6 M and 
- 1.1 x 10 6 M Q ) at a redshift of just 1.25. What is 
disturbing is the own existence of this bound and the 
related unsatisfactory behavior of the merger rate at 




Fig. 1. — Instantaneous destruction (or true merger) rate at the present time for an object with mass M/M Q equal to 10 8 
(solid lines), 10 10 (dashed lines), 10 12 (long-dashed lines), and 10 14 (dotted-dashed lines) as a function of the mass AM of 
the merger partner for the CDM (a) and n = — 2 power-law (b) power spectra, both normalized to a(8h^ 1 Mpc) = 0.67 
in an Einstein-de Sitter universe with h = 0.5. In thin lines the merger rates drawn from the LCa model. 



moderate AM/M. These shortcomings do not seem 
to be inherent to the basic assumption of the model, 
that is, the possibility to effectively distinguish be- 
tween accretion and merger. They could be partly 
due to the fact that the peak model is only a good 
approximation for the gravitational growth of mas- 
sive objects. (At any given 5 C , the larger the scale, 
the higher the peak amplitude as compared to the 
typical fluctuation, ao(R)- And since the higher the 
peak, the more spherical its shape and the more negli- 
gible the shear caused by the surrounding matter, the 
better its gravitational collapse is described by the 
spherical approximation.) But what the most likely 
causes them is the approximation used in equation 
(1). Indeed, some preliminary calculations seem to 
indicate that a more accurate approximation for the 
nesting correction should notably improve the perfor- 
mances of the model. In the following, we will not 
insist anymore on the limited domain of consistency 
of the present model (in its current version) and focus 
on the analysis of its predictions for large enough M. 

From the confluent system formalism making the 
natural distinction between accretion and merger we 
can also calculate the transition rate including both 
kinds of captures, mergers and accretion. Equation 



(5) gives the rate at which objects of mass M merge 
and are destroyed giving rise to objects of mass M' . 
But from the viewpoint of these latter objects the 
event can be seen either as accretion, if they can be 
identified with the initial most massive partner, or 
as a true merger, if they cannot. Thus, the instan- 
taneous accretion+merger, or simply capture, rate 
for final objects of mass M' per specific infinitesimal 
range of the captured mass M (M < M') is 



r c (M' <- M,t) 



r d (M -> M',t)N(M,t) 
N(M',t) 



(6) 



with r d given by equation (5). This capture rate is not 
directly comparable to the LCa merger rate plotted in 
Figure 1 since it refers to the mass M' of the final ob- 
ject instead of to the mass M of the initial object. But 
it can be readily compared to the capture rate drawn 
from the LCa merger rate through the same relation 
(6). This comparison is shown in Figure 2. Once 
again we find a markedly different behavior between 
the solution of the two models. This difference is 
also due to the distinct viewpoint adopted in the two 
models. The capture rate obtained from our model 
making the distinction between accretion and merger 
only includes, in the case of accretion, those objects 




Fig. 2. — Instantaneous capture (or accretion+ merger) rate at the present time for a final object with mass M /Mq equal 
to 10 14 (the lowest curve on the left), 3 X 10 13 , 9 X 10 14 , 2.7 X 10 15 , and 8.1 X 10 15 (the highest curve on the left) as a 
function of the mass M of the captured objects for the same cosmogonies (a) and (b) as in Fig. 1. In thin lines the capture 
rates inferred from the LCa merger rates. 



which are accreted, obviously not the accreting ones 
identified to the final objects of mass M', and hence, 
not considered as captured by them. In contrast, all 
objects are counted in the LCa capture rate (in par- 
ticular, those seen as accreted as well as those seen as 
accreting from our point of view) because there is no 
identification in the LCa model between the final and 
any of the initial objects partaking of any capture. 
This clearly explains the difference between the two 
capture rates at captured masses M close to M'. On 
the other hand, the two models should show a similar 
behavior at small M, in full in agreement with what 
is found. 

The previous destruction and capture rates are per 
specific infinitesimal range of mass of the final or cap- 
tured objects, respectively. To derive the respective 
global rates we must simply integrate equation (5) 
over M' 

/>oo 

r d (M,t)= r d {M -» M',t)dM', (7) 
Jm 

and equation (6) over M 
l~M' 

r c (M\t)= r c (M' <- M,t)dM. (8) 
Jo 



Note that these global rates are obviously positive 
down to much smaller masses M or M' than the re- 
spective specific ones for the whole range of the ex- 
tra variable. It is worthwhile mentioning that these 
global merger and capture rates cannot be compared 
with the respective quantities in the LCa model be- 
cause these cannot be calculated in that model (the 
corresponding integrals diverge). This is an impor- 
tant drawbak of the PS approach related to the fact 
no distinction is made between merger and accretion. 

Let us now turn to the formation rate. The den- 
sity of objects that form (as the result of mergers) in 
the interval of time dt is traced by the net density 
of non-nested appearing peaks yielded in the corre- 
sponding range —dS. More exactly, we are interested 
in calculating the net density of non-nested peaks 
with S c on scales between R and R + dR that just ap- 
pear (i.e., fusion minus fission) in the preceding range 
—dd, Nf(R,5 c )dRd5. This net density is derived 
in Appendix C. Now, by dividing it by N(R, S c ) dR 
we are led to the net conditional probability that a 
non-nested peak with scale R appears between 5 and 
S c — dS. Therefore, the instantaneous formation rate 




at t of objects of mass M is 



d5 r 



dt 



(9) 



with R and 5 C on the right hand side written in terms 
of M and t, respectively, through equations (3) and 
(4). This formation rate is plotted in Figure 3. No 
similar quantity is provided by the LCa model. 

Finally, by multiplying equation (5) by AM — 
M' — M we obtain the specific rate at which the mass 
M of objects is increased at the time t owing to (true) 
mergers and, by integrating this latter function over 
M', the instantaneous typical mass increase rate ow- 
ing to mergers for objects of mass M 



/•OO 

(M,t)= AM r d (M —* M', t) dM' . 
Jm 



(10) 



Likewise, from equation (6) we can obtain the instan- 
taneous typical mass increase rate for objects of final 
mass M' owing to accretion and merger, 



,(M',t)= / 
Jo 



M r c (M' <— M, t) dM. (11) 



As shown in Appendix D, the idea of a transition 
rate, by accretion, similar to those given by equations 



(5) or (6) is meaningless because accretion is a con- 
tinuous instead of a discrete process of mass increase. 
In contrast, the mass increase rate by accretion does 
make sense. The instantaneous mass accretion rate 
for objects of mass M follows from the instantaneous 
scale increase rate of the corresponding peaks as they 
evolve along continuous and derivable trajectories in 
the S vs. R diagram. This latter rate depends on the 
particular value of the scaled Laplacian x of the peak 
which is being followed. However, since we are inter- 
ested in the typical mass accretion rate for objects of 
mass M disregarding any other property we must av- 
erage over x. This is done in Appendix D. The result 
is 



,(M,t) 



1 



dM 



dS c 



dt 



(12) 



(x)a 2 R dR 

with R and S c on the right hand side written in terms 
of M and t, respectively, through equations (3) and 
(4) and (x) the average mentioned after equation (1). 

In Figure 4, we plot the total mass increase rate of 
objects of mass M at t, r mass (M,t) = r^ ass (M 7 t) + 
r mass (-^i t) j an d compare it with that resulting from 
the LCa model, i.e., the integral over AM of their spe- 
cific merger rate times the corresponding increase in 
mass. (Note that this latter integral does not diverge 
contrary to those involved in the calculation of the 
total merger and total capture rates.) As expected, 




Fig. 4. — Instantaneous total mass increase rate for the same epochs and masses as in Fig. 3. The contributions from 
mergers and accretion (dominating at the lower and higher mass ends, respectively) are represented in dashed lines. In thin 
lines the total mass increase rates inferred from LCa merger rates. 



the two models show, in this case, similar qualita- 
tive behaviors. They only differ quantitatively in the 
regime dominated by mergers. This departure is just 
due to the differences in the two merger rates at large 
AM/M (see Fig.l). 

4. GROWTH CHARACTERISTIC TIMES 

From the meaning of the global merger rate, equa- 
tion (7), we have that the density N sur (t) dM of ob- 
jects surviving (i.e., having not merged but just ac- 
creted) until the time t from a typical population with 
masses in the range between Mq and Mq + dM at 
to < t is given by the solution of the differential equa- 
tion 

^L = -r d [M(t),t]N sur (t), (13) 

with initial condition N sur (t a ) = N(M ,to)- In equa- 
tion (13) and hereafter, the function M(t) is the typi- 
cal mass at t of such accreting objects, approximately 
given by the solution of the differential equation 

^=r a mass [M(t),t], (14) 

with initial condition M(to) — Mq. Indeed, from 
equation (12) we have that the typical mass of ob- 
jects evolving by continuous accretion increases with 
time according to equation (14). This expression is 



only approximated because the average over x leading 
to equation (12) presumes all objects with the same 
mass, while the mass M(t) is just the so estimated 
typical mass of objects at t. Hence, the larger the in- 
terval of time spent since to when objects have really 
identical mass the poorer will be the approximation. 

The solution of equation (13) is 

N sur (t) -7V(M ,to)exp|-^ r d [M(t'),t']dt'^. 

(15) 

Hence, by defining the typical survival time t sur of 
objects with masses between Mo and Mo + dM at to 
as the interval of time since to after which their initial 
density is reduced (owing to mergers) by a factor e, 
we are led to the equality t sur = td — to, with the 
destruction time td(M ,to) given by the solution of 
the implicit equation 

rt d (M ,t ) 

1 = / r d [M{t'),t']dt' . (16) 

J t 

In addition, the typical mass accreted by those objects 
until they merge and disappear is M[td(Mo, to)] — Mo- 
(Caution: what LCa called survival time rather corre- 
sponds to what here is called destruction time; here- 
after we assume these authors using the same notation 
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as ours.) 

In a fully similar manner we can infer the typical 
age of objects with masses in the range between Mo 
and Mo+dM at to, that is, the typical interval of time 
since the last merger giving them rise. The density 
N pre (t)dM of these objects pre-existing (i.e., having 
just accreted matter since then) at a time t < to is 
given by the solution of the differential equation 

= rf[M{t),t]N[M{t),t] - r d [M(t),t]N pre (t) 

(17) 

with initial condition N pre (t ) = N(M ,to). The so- 
lution of equation (17) is 

N pre (t) =N(M Q ,to)K- 1 (t,t ) 

1 -/ v liM lr '[ M<, '»^'i A ' (I ' w } (18f '» 

K(t, to) = expj- J r'lJf ({),£] dA (18b) 

Thus, by defining the typical age t age (M ,to) of 
objects with masses between M and M +dM at to as 
the interval of time until to before which their density 
(owing to their progressive formation) was a factor e 
smaller, we arc led to the equality t age =to~tf, with 
the formation time i/(Mo, to) given by the solution of 
the implicit equation 



N pre [t f (M ,t )} 
N(M ,t ) 



c. 



(19) 



And the typical mass accreted by these objects since 
they formed is M - M[t f (M , t )]. 

Finally, from the typical age and surviving time 
of objects with masses between M and M + dM 
at to calculated above, we can readily calculate their 
typical lifetime or intermerger period. This is simply 



tuf e (M ,to) = t age (Mo,to) +t sur (M ,to) 



(20) 



= t m (Mo,to)-t f (M ,to). 

And connected with this latter quantity, there is 
the total mass typically accreted by those objects 
during their whole life, given by M[t m (Mo,to)] — 
M[tf(M ,t )]. 

In Figures 5 and 6 we plot the typical age and 
survival time, respectively, of objects at two different 
epochs as a function of their mass M. For comparison 
we also plot the corresponding time estimates drawn 



from LCa model. As mentioned, our natural defini- 
tions of the typical age and survival time (i.e., the 
time spent since the previous merger when the object 
was formed and the time spent until the next merger 
when the object will be destroyed, respectively) are 
not possible in the PS approach not making the dis- 
tinction between merger and accretion. What LCa 
call the age and survival time of an object of mass 
M are the time spent since the mass of some parent 
object reached half that mass and the time required 
by the object to be found inside an object of double 
mass, respectively. (LCa give the time distributions; 
what are plotted in Figs. 5 and 6 are the correspond- 
ing median values.) 

Our age shows the same slightly decreasing trend 
with increasing M and the same strongly increasing 
trend with increasing time as LCa's. The only dif- 
ference is that our age is slightly smaller than LCa's. 
(This trend is reversed at very large masses, for a 
given redshift, but this is likely due to the poor fit 
to the PS mass function yielded there by our mass 
function. Note that a mass of 10 12 M Q at z = 2 
is comparable to a present mass of two or three or- 
ders of magnitude higher for the n = — 2 power-law 
and the CDM spectra, respectively.) In contrast, our 
survival time show a very different, in fact opposite, 
behavior from LCa's. Our result that the larger the 
mass of the object, the larger its survival time is just 
what one would expect in hierarchical clustering dis- 
tinguishing between accretion and merger. Indeed, at 
a fixed time, massive objects are rare and they hardly 
merge, while light objects are more frequent and can 
more easily merge (see Fig. 3). Moreover, as the time 
increases objects of a given fixed mass become more 
and more frequent, hence, their survival time dimin- 
ishes as also found. 

To better understand these comparative behaviors 
we also plot, in Figures 5 and 6, two new times, 
hereafter called half-mass-accretion time and double- 
mass-accretion time, defined a la LCa but in the con- 
text of our model. These are the interval of time spent 
since the mass of an object was typically half its cur- 
rent value and the interval of time required by an 
object to typically double its mass, respectively. As 
these new times refer to the typical mass evolution 
of given objects the evolutionary process involved can 
only be accretion. So they are given by the function 
M(t) solution of equation (14). It is worthwhile not- 
ing that these new times are, in some sense, orthogo- 
nal to our age and survival times because the former 
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Fig. 5. — Typical age (solid lines) and half- mass-accretion time (dashed lines) of objects of mass M in units of Mq = 
10 12 Mq for the quoted values of the redshift corresponding to the same cosmogonies (a) and (b) as in the preceding figures. 
LCa's median age is plotted for comparison (thin dotted lines). 




M/M r , M/M r , 

Fig. 6. — Typical survival time (solid lines) and double-mass-accretion time (dashed lines) of objects of mass M in units of 
Mq = 10 12 Mq for the quoted values of the redshift correponding to the same cosmogonies (a) and (b) as in Fig. 4. LCa's 
median survival time is plotted for comparison (thin dotted lines). 
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only depend on accretion while the latter mostly do 
on mergers. 

From Figure 5 we see that objects found at a given 
rcdshift typically form at a lookback time of the order 
of (just slightly smaller than) that required for them 
to accrete half their mass. (Note that our age takes 
into account that many similar objects that formed 
earlier or later than this epoch have not survived un- 
til the redshift considered.) So mergers play an im- 
portant role in the past evolution of these objects, 
but the mass increase experienced by them from the 
time they typically formed is close to a factor two. 
It is therefore not surprising that the age estimate a 
la LCa (accounting for both accretion and mergers) 
is also close to that time. Notice, however, that the 
similarity between the half-mass-accretion time and 
the age as defined in the present work is fortuitous 
since both time estimates are rather orthogonal. So 
is the similarity between our age and LCa's. 

From Figure 6 we see that the time required by 
very massive objects at a given epoch to be destroyed 
is much greater than that required for them to double 
their mass by accretion, while the opposite is true for 
very light objects. In other words, there is a gradation 
from frequent light objects where the mass increase is 
dominated by mergers to rare massive objects where 
it is by accretion (see Fig. 3). It is therefore not 
surprising that the double-mass-accretion time (only 
accounting for accretion) and the LCa survival time 
(accounting for accretion and mergers) notably devi- 
ate from each other for very light objects while they 
yield similar results for very massive ones. Note that 
we should not expect a very good agreement there, 
either, because the larger the resulting value of the 
time estimate, the more marked the effects of: 1) the 
approximation used (eq. [14]) to calculate accretion- 
times and 2) the difference between the two couples of 
time estimates; ours reports on "the time required for 
objects to change their typical mass by a given fac- 
tor", i.e., they involve averages along R at fixed <Ts, 
while LCa's reports on "the typical time required for 
objects to change their mass by the same given fac- 
tor", i.e., they involve averages along 5 at fixed R's. 
In any event, the opposite behavior shown by our sur- 
vival time as compared to LCa's is also clarified. 

5. DISCUSSION AND CONCLUSIONS 

The idea that the present clustering model relying 
on the peak model ansatz can provide a good descrip- 
tion of hierachical clustering is at variance with the 



rather extended opinion that peaks are not good seeds 
of bound objects. Detailed TV-body simulations of the 
gravitational evolution of individual density maxima 
seem to show, indeed, that if their initial amplitude is 
small they can be disrupted before collapsing (van de 
Weygaert & Babul 1993). While Katz, Quinn, & Gelb 
(1993) have found that even high amplitude peaks 
seem to be poor tracers of bound objects, individu- 
ally as well as statistically. The reason why this would 
be so is that peaks are not spherical in general, nor 
are they isolated. So their non-linear evolution will 
markedly deviate from the spherical approximation at 
the base of the peak model ansatz. However, any con- 
clusion drawn from the previous TV-body simulations 
dealing with peak statistics is somewhat precipitated 
because of the unknown effects of the use of unap- 
propriate filtering windows and/or R(M) and S c (t) 
relations, and the lack of any correction for the cloud- 
in-cloud effect. In fact, by using the accurate peak- 
patch formalism, Bond & Myers (1996) have found 
that there is actually a good correspondence between 
peaks and high mass objects. This is well understood: 
the higher the peak relative to the typical density fluc- 
tuation on that scale, the more spherical, and the less 
important the shear caused by the surrounding den- 
sity fluctuations, particularly if the power spectrum is 
steep enough (<r (M) flat enough) to guarantee that 
nearby density fluctuations on similar scales are typi- 
cally negligible. This has also been recently confirmed 
by Bernardeau (1994) who has rigurously shown that 
the evolution of high amplitude peaks (with v ;> 2) in 
a Gaussian random density field is correctly described 
by the spherical model provided only the logarithmic 
slope of the power spectrum is smaller than n = — 1. 

A fundamental result of paper I was that the self- 
consistency of the peak model ansatz, in an extended 
version presuming some effective distinction between 
accretion and merger, completely fixes the filtering 
window and the R(M) and S c (t) relations except for 
a couple of free parameters governing the dynamics of 
collapse. After calculating the density of peaks at a 
fixed S c and performing the appropriate correction for 
the cloud-in-cloud effect we obtained a mass function 
which is in overall agreement with iV-body simula- 
tions (since in agreement with the PS mass function) 
for appropriate values of those parameters. 

In the present paper we have calculated the growth 
rates and characteristic times of relaxed objects pre- 
dicted by the model. As far as we can tell from 
the comparison with the LCa model not making the 
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distinction between merger and accretion, our model 
yields a reasonably good detailed description of hi- 
erarchical clustering for massive objects. Its unsat- 
isfactory behavior at small M is likely due to one 
practical although rather crude approximation used 
in paper I to correct for the nesting among peaks at 
a fixed density contrast. We are currently working on 
an improved version of the model based on a more 
accurate correction. 

Let us finally insist on the great potential of the 
approach followed here. The violent relaxation that 
accompanies the sudden large mass increase of a re- 
laxed system is expected to leave perdurable imprints 
in its morphological properties as compared to the ef- 
fects of a quasi-continuous one. A detailed clustering 
model should therefore make that important distinc- 
tion. This can be readily carried out in the extreme 
cases of very large or very small mass captures, but 
in the intermediate regime it is more problematic. Do 
we have we the right to assume some definite effec- 
tive division between tiny and notable captures? Is it 
unique and where is it? The assumption that there 
is, indeed, a definite division between accretion and 
merger yields, in the well motivated framework of the 
peak theory, one unique, self-consistent, reasonably 
good clustering model. Thus, the possibility of car- 
rying out such a division in a self-consistent manner 
seems proved. The comparison with iV-body simu- 
lations should tell us where the effective division be- 
tween merger and accretion is, opening in this way the 
possibility of better describing the growth of cosmic 
objects and of understanding their observed proper- 
ties. 

We are grateful to the referee, Cedric Lacey, for 
his fruitful criticisms. This work has been supported 
by the Direction General de Investigation Cientffica 
y Tecnica under contract PB93-0821-C02-01. 
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APPENDIX A: PREPARATORY RESULTS 

A.l. DISAPPEARING PEAKS 

From the Taylor series expansion of the gradient of the density contrast smoothed on scale R, rj, around the 
location r p of a neighboring peak we have (to first order in |r — r p |) 

Vi ~ (r ~ r p )i, (Al) 

with Aj > 0, the eigenvalues of the second order cartesian derivative tensor Qj changed of sign evaluated at r p or 
r. Since there is at most one peak in the neighborhood of any point the density of peaks around a point r is 

6 (3 \r-r p ) = |A 1 A 2 A 3 | ( 5 ( S- (A2) 

But r]i and Aj are random Gaussian variables. So the typical density of peaks on scale R around an arbitrary point 
r is given by the mean 

(<5( 3 )(r-r p )) HIAxA^l^r,) (A3) 

for the joint probability of the random variables involved, with Aj strictly positive, evaluated at the arbitrary point. 
This is the scheme followed by BBKS for obtaining the density Af p k(v, R) dv of peaks with scaled density contrast 
v in an infinitesimal range. 

According to the identification criterion among peaks on contiguous scales, a peak at r p on scale R will disappear 
in reaching the scale R + AR, with AR positive and arbitrarily small, provided only this is the first scale larger 
than R with no peak in the neighborhood |r — r p | < O(AR) of the former peak. From the Taylor series expansion 
of rj at the nearest point r p > with r] = on scale R + AR around the peak r p on scale R we have 

(V-r p )i«^Ai2. (A4) 

Thus, on the new scale R + AR there is some peak in the neighborhood |r — r p \ < O(AR) of the old peak on scale 
R provided only that all Aj are of order unity (this also guarantees that all Aj are positive at r p > on scale R + AR). 
Strictly, some A, could be of order AR or smaller, if the corresponding darji were too. But the probability for 
this to happen is negligible as compared to the more general preceding case. Therefore, for that condition to be 
broken for the first time at R + AR, some eigenvector Aj must become of order AR or, equivalently, must vanish 
at R ± AR. 

Consequently, the density of peaks on scale R at an arbitrary point r disappearing at R' in the neighborhood of 
R, with R' > R, is given by the mean {5^{r — v p )), with all eigenvalues Aj positive, as in BBKS, and the smallest 
one satisfying the relation 

R '- Ra wk\- (A6) 

But condition (A5) implies, at the same time, that the smallest eigenvalue Aj also vanishes in the neighborhood 
of r p on scale R. Thus, in neglecting second order terms in equation (Al) the term in that component Aj must 
be neglected, too. Accordingly, were we interested in calculating the density of peaks disappearing at R' in the 
neighborhood of R, we should take for that component of r — r p , instead of (Al), the relation 

dnVi - d R f]i{v p ) w -d R \i (r - r p )j, (A6) 

with duT]i(r p ) evaluated at the peak on scale R and the remaining variables evaluated at the arbitrary point r. 

Following these schemes it can be shown (Manrique 1995) that the density of peaks at R with v in an infinitesimal 
range and all eigenvalues Aj of zeroth order in AR is equal to the total density Af p k(v, R) dv calculated by BBKS 
to zeroth order in AR, while the density of peaks with some Aj of first order in AR contribute to that total density 
with first order terms in AR. Therefore, the density of peaks at 5 with R in an infinitesimal range, N p k(R,S) dR, 
calculated in paper I from the density Af p k {v, R) dv is indistinguishable to leading order in AR (or in AS for each 



Growth History of Objects 



17 



given trajectory S(R)) from the density of peaks at 5 robust enough for their continuity to be guaranteed in some 
neighborhood AR (or AS). 



A. 2. DENSITIES OF PEAKS WITH SPECIFIC VALUES OF RANDOM VARIABLES 

In next Appendixes we will be concerned with densities of peaks with values of the random Gaussian variables 
with null mean vq = v, v\ = x = — (Ra 2 ) d^aov) (this is the same variable x defined in BBKS; see paper I), and 

d R [o2{i-i)K- 2 Vi-i] (i>2) (A7) 



a 2t R 1 - 1 

on scale R in infinitesimal ranges. All these densities can be inferred from the density of peaks with v and x on 
scale R in infinitesimal ranges (BBKS) 



exp(— v 2 12) 

Af pk (v,x,R)dvdx = {2n)2m [2?r(1 _ 72)]1/2 f(x) exp 



(x - ~fv) 2 
"2(l-7 2 ) 



dv dx (A8) 



with f(x) a function given by BBKS (eq. [A15]), 7 = a\ / {aoa?) , and R„ = \Z3<r \j 02- One must simply apply the 
recursive relation 

Af P k{vo,vi, v n , R) dv dvi ...dv n = N p k{vo,v\, v n -i,R) dv Q dvi, ...dv n -iV(v n ,R\v , vi, ...,v n -i,R) dv n , (A9) 

with V(w, R\vo, v\, ...,v n -i, R) dw the conditional probability of finding the value of w in an infinitesimal range 
given that vo,v\,..., v n -\ take some given values, 



V(w,R\v ,Vi, ...,v n -i, R)dw = —=L exp 

V2na w 



(w — w) 2 



, , dw (AlOa) 

w = M„M„"JV T ^=<to 2 >-M„M„ 1 M^, (AlOb) 

where M Kt and M„„ are the 1 x n and the nx n correlation matrixes of w with Uj, (w Vi), and Vi with themselves, 
(vi Vj), respectively, V stands for the 1 x n matrix of components Vi, and T denotes transpose. Note that, although 
equation (A9) involves, in principle, the conditional probability for peaks, the conditional probability (A10) is for 
simple points. The reason for this is that the condition for a point to be a peak refers to the values of variables 
r)i, x, y, and z (see BBKS for the definition of y and z). While none of the variables vi defined above correlates 
with any of these variables except for x. Consequently, the conditional probability of finding any specific value of 
Vi does not depend on whether we are dealing with a peak or simply a point with the same given (positive) value 
of x. 

The correlations among variables Vi take, for the Gaussian window, the general form 

2 

( Vi Vj) = (-iy+3+^+So 3 _^±J_ ( AU ) 

{hi > 0), with Sij the Kronecker delta. We must remark that, although the notation used in the preceding 
equations presumes all variables defined on scale R, they are also valid for variables defined on different scales by 
just appropriately changing the values of the spectral moments involved. In particular, were any variable, say Vi, 
defined on another scale R' , one would be led to just the same expression for the correlations as in equation (All) 
but with cjj and Ui + j replaced by a\ = <Ji(R') and a i+ j h = <Tj + j{[0.5 (R 2 + R' 2 )] 1 / 2 }, respectively. 

Thus, from equations (A8) to (All) one can readily calculate the density function Af p k of peaks at a fixed scale 
R per infinitesimal ranges of v, x, and any other set of the previous variables Vi (i > 2). And by integrating it over 
variable v with the constraint (see paper I) 

— <v< Vx — R AR, (A12) 

(T (Tq (To 
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we can infer the density function N pk of peaks at a fixed 8 per infinitesimal ranges of R, x, and the same set of 
variables Vi (i > 2). Consequently, such density functions N pk are just equal to the corresponding ones M pk times 
x a 2 R/a . 

Moreover, following the same scheme starting from the conditional density analogous to the normal density in 
equation (A8) but for peaks subject to any given constraint one is led to the same relation as above between couples 
of conditional density functions N pk and N pk per infinitesimal ranges of any set of variables Vi (i > 2) subject to 
that constraint. 



A. 3. NESTING PROBABILITIES 

Collapsing clouds associated with non-nested peaks with fixed 8 yield a partition of space which makes the mass 
function of objects at t be correctly normalized (see paper I). This implies that the volume fraction occupied by 
disjoint backgrounds with 8 on filtering scales between R' and R' + dR' or, equivalently, the probability to find a 
point in any such backgrounds is M(R') p~ x N(R', 8) dR' . Therefore, 

P ,M W )*=B»( W )*f*a. ( Ai3) 

P N pk (R,d) 

gives the approximate probability that a typical (non-fissionable) peak with 8 on scale R is nested within some 
non-nested peak with identical density contrast but on a scale between R' and R' + dR' (R < R'), hereafter simply 
called the (differential) nesting probability of a peak. (Remember that the density of peaks that fission in the next 
dR or — dS is a higher order differential.) This was used, in paper I, to derive equation (1). Likewise, the nesting 
probability of peaks with given values of variables 8, x, and any set of variables Vi (i — 2, n with arbitrary n) is 

p^r'aira W»' M (^) auw *wp/ N pk (R,x,v 2 ,...,v n ,S\R',S) 

P(R ,8\R,8,x,v 2 ,...v n )dR = N(R,S)dR F — , (A14) 

p N pk (R,x,v 2 ,...,v n ,8) 

with N pk (R, x, v 2 , v n , S\R', 8) the conditional density function analogous to N pk (R, 8\R', 8) in equation (A13) but 
per infinitesimal ranges of the extra variables x, v 2 , V3, v n . The last factor on the right-hand member of equation 
(A14) satisfies the relation (see Appendix A. 2) 

N pk (R,x,v 2 ,...,v n ,5\R',S) _ 7V, R'\u, x, v 2l v n , R)N pk {v, x, v 2 , ...,v n ,R) x^R _ V(v' 7 R'\v,x,v 2 , ...,v n , R) 
N pk {R,x,v 2 ,...,v n ,8) ~ V(v',R')N pk (v,x,v 2 ,...,v n ,R)x^R ~ V{v',R>) 

(A15) 

with v' = S/a' , V{y' , R'\v, x, v 2 , v n , R) dv' defined in equation (A10), and V{v,R)dv the Gaussian probability 
to find the scaled density contrast on scale R between v and v + dv. Notice that, for identical reasons as for the 
conditional probability in equation (A10), the conditional probability given in equation (A14) applies, in fact, to 
points. Note also that we are using the same convention for the notation of the conditional probabilities as for the 
density functions: they are denoted by a caligraphic capital p (in contrast with the notation used in BBKS) when 
they are per infinitesimal range of v at a fixed R, and by a roman capital p when they are per infinitesimal range 
of R at a fixed 8. Finally, it is worthwhile mentioning that, contrary to what is suggested by the present notation, 
we can write v or 8, indistinctly, when specifying the condition. 

From equation (A10), it is clear that the nesting probability given by equations (A14) and (A15) will depend 
on variables Vi provided only that these variables have non-null correlations with v (or 8) and x on scale R or 
any explicit variable correlating with them. This is what happens with the variables Vi defined in equation (A7) 
involving the different order scale derivatives of the density contrast. (This is readily understood from the Taylor 
series expansion of v on scale R! around R: the probability that a point has v on scale R! depends on the value 
of all scale derivatives of the density field on scale R.) Thus, to accurately infer the density of nested peaks 
in any given peak population one must calculate the distribution of these infinite variables in that population. 
The only noticeable exception concerns the typical population of peaks. The density of typical peaks at 8 with 
values of Vi (i = 1, n) in infinitesimal ranges, N pk (R, x, v 2 , ...v n ,8) dRdxdv 2 ...dv n , times the nesting probability 
P(R' ,S\R,8 7 x,v 2} ...v n ) dR' integrated over any subset of variables Vi coincides with the product of these two 
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functions without the explicit dependence on the integrated variables, as readily seen from equation (A14) valid 
for any arbitrary values of the subindexes. In particular, the integral over all Vi for i > 1 is equal to the product of 
N pk (R,S)dR times the reduced nesting probability P(R' ,5\R,6) dR' given in equation (A13), which justifies this 
latter expression and equation (1). 

The conditional probability P(v', R'\v, x, v 2} V3, R) dv' extended to the infinite set of variables Vi defined in 
equation (A7) can be obtained according to equation (A10). After some lenghty algebra using intermediate variables 
which only correlate with themselves (found by means of the Gramm-Schmidt method) we arrive to the expression 



V(v' ', R'\v, x, V2,v 3 , R) dv' 



1 



2ir ov< 



exp 



dv 1 



(A16) 



with 



v = a is + ol\X + a 2 v 2 + 



OLi = ^2 B 3 & J-* X! ^ k 3- k < V Q V k> U l' = 1 ~ ^ B l X! ^ k 3~ k < V V k> (A17) 

j=0 fc=0 j=0 L k=0 

(i > 0), where a prime denotes the scale R', the correlations in angular brackets are given by equation (All), and 
cocffccients Bj and (5 k j- k are defined as 



j-i / i-i-fc 

(Bj)' 1 = 1 + (5 k j_ k I (5 k j-k + 2 fa- 1 1 < Vk v i~ l > 



Pij-i — ^ C, 



» 3~* 



< Vj v k >, 



fc=0 x Z=0 7 fe=0 

(j > 1 and i 7^ j) in addition to Bq = 1 and /3 i0 = 1- In equations (A18) we have used the notation 



Cfl = B i Pk i-k 



3-1 

^ij-i = Bi(5 k (5i , 
l=k 



(A18) 



(A19) 



for i < k and C k j _ i = C l k -_ k for % > k. 



From the general expressions of w and a^, equation (A10), it can be shown, through the intermediate use of 
variables defined as v but for just the first arbitrary n + 1 variables Vi, that 



( VVl ) = (v> Vl ) = (_i)i+i+«oi 

C CT2j 



(A20) 



(i > 0) which implies, on its turn, a similar relation for any (finite or infinite) linear combination of Vi, in particular, 

(V 2 ) = (v'y) = l-a 2 v ,. (A21) 

On the other hand, taking into account that the correlation between two variables is equal to the integral of the 
product of their Fourier transforms taken at r = and defining the new variable x, as — (cr 2 R)" 1 dji(<j' v) we have 



(xvi) = 

C2 



f a 2i+l 1 o \ ,_ \ . / , n«5 i0 °0 °2(i+l) ,_ 
- — p^fl Wl/Vj) + (-1) (uv i+1 



(i > 0) and 



(x v) 



By substituting the correlations in equation (A20) into equation (A22) we obtain 

(xvi) = 



(A22) 



(A23) 



(A24) 
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(i > 0) which also implies 

(xu)=0 (x 2 } = 0. (A25) 

Thus, (P 2 ) and ov do not depend on R (see eqs. [A21] and [A23]) which ultimately implies that a v t is null and 
{p ) equal to unity. Indeed, according to its definition, equation (A17), ov is null for R' = R and since it does not 
depend on R it is necessarily null for any value of R'. Then, equation (A16) leads to 

r(u',R'\u,x,v 2 ,v 3 ,...,R)du' = 8{v' -v)dv'. (A26) 



This result is not surprising since fixing the values of the density contrast and every order scale derivative of it on 
a given scale R automatically fixes, through the Taylor series expansion of 8 as a function of the filtering scale, 
the value of the density contrast in any other scale R! . Finally, by substituting V given by equation (A26) into 
equation (A15) and the latter into (A14) we arrive to the following expression for the nesting probability 

P(R', S\R, 8, x, v 2 ,v 3 , ...) dR' = P(R', 6\v) dR' = g) dR' , (A27) 

with v' = 5/a' . 



APPENDIX B: NET DENSITY OF PEAKS BECOMING NESTED 
The density of peaks at Sf = 8 — Ad, with AS positive and arbitrarily small, on scales between R f and Rf + dR f 
and variables Xf,v 2 f,v 3 f,... in infinitesimal ranges, which result by continuous evolution (hence, without fissioning) 
from peaks at S on scales between R and R + dR and x,v 2 ,v 3 , ... in infinitesimal ranges is 



with 



N*%(Rf,Xf,v 2 f,v 3 f, ...,Sf) dRfdx f dv 2f ... = N pk (R,x,v 2 ,v 3 , ...,S) dRdxdv 2 dv 3 

AS 

I 

AS 



R f w R + 



xa 2 R 



v 2 



04 

e>2 



V 2 f ~ v 2 



&6 



V3 h v 2 \ — - 

(T 4 



AS 
a 2 J xa 2 

1 



R 2 



— + (E 

xa 2 \ ^ 



xa 2 R 



div 2 



AS 
xa 2 R 



(Bl) 
(B2) 



(B3) 



to first order in AS. Equation (Bl) states that the density N p k{R 7 x, v 2 , v 3 , 8) dRdx dv 2 dv 3 ... of peaks is 
conserved through continuous evolution from 8 to Sf. Equations (B2) and (B3) give the shift from 8 to Sf in 
the values of all the relevant variables. It is important to outline that we need to know, indeed, the values of 
all variables R and Vi (i > 1) at 8f in order to calculate the density of evolved peaks which are nested since, as 
explained in Appendix A. 3, the nesting probability depends explicitly on all these variables differently distributed 
in than in N pk at Sf. Equation (B2) arises from the derivative dR/dS along continuous peak trajectories, equal 
to — (xa 2 i?) _1 (see paper I). While the shift in the variables Vi, equations (B3), is equal to the sum of two terms: 
one coming from the scale derivative of each particular variable, and a second one coming from the scalar product 
of its spatial gradient times the shift in position of the new peak relative to the old one (eq. [A4]). However, the 
nesting probability does not depend on the variables dsJIi since these variables do not correlate with v, x, v 2 , w 3 , ... 
(see Appendix A. 3). Thus, in averaging below over all variables, these second terms will contribute with a null 
mean (the distribution of dsJli for peaks is the same as for arbitrary points). Consequently, we can drop these 
second terms which is equivalent to taking an effective location for each evolved peak equal to the mean expected 
value, that is, the same location as the original peak. 
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From equation (Bl) we have that the density of peaks at 5 per infinitesimal ranges of R and x, V2, V3,... which, 
after evolving to Sf, are found to be nested (although not necessarily become nested) into non-nested peaks with 
scales between R' and R' + dR' (R < R f < R') is 

N^ st (R -> R' , x, v 2 , v 3 , . . . , 5 -> S f ) dR dR' dx dv 2 dv 3 . . . 

= N p k(R, x,V2,V3, ...,5) dRdxdv2 dv 3 ... P(R! , 5f\Rf, Sf,Xf,V2f,v 3 f,...) dR' , (B4) 

with Rf, Xf, V2f, vsf,... on the right hand side in terms of R, x, v 2l v 3 ,... and AS through equations (B2) and 
(B3), and the specific nesting probability P given in Appendix A. 3. To obtain N^ st (R — > R',S — > Sf)dRdR' 
giving the density of peaks at 5 per infinitesimal range of R which after evolving into Sf are found to be nested into 
non-nested peaks on scales between R' and R' + dR', we must integrate equation (B4) over variables x,v 2 ,v 3 , .... 
Given the simple expression of the nesting probability appearing in equation (B4) in terms of the variable v (eq. 
[A27]) it is convenient to first transform v 2 and v 3 to v and x. This can be done by repeated application of the 
scheme given in Appendix A. 2. Taking into account the correlations (A20) and (A21) and the fact that cr 2 , = 
(see the discussion after eq. [A23]) we first obtain 



e 



N pk {R,x,v,R' ,S) dRdxdv ~ N pk {R,x, 8) dRdx = dv, (B5) 

V27T(7p 



(B6) 



with N p k(R, x, S) = N p k{v, x, R) x a 2 R/<tq in terms of the density function given in equation (A8) and 



S 2,2 

ft =a P0 Yolv\x a 9 = 1 - e 



(I-7 2 n) 2 



1 — 7 



— +r\ 



era 

e(l-7 2 ri) 67(1 -n) 

apo = : 9 ota 1 — ^ — 

1—7^ 1—7^ 

with e = al h /(ao<r' ) and n = o'ih a o/( <T oh cr i) already used in BBKS. Then, from correlations (A24) and (A25) 
we are led to 

N pk (R, x, v, x, R! , S) dRdxdv dx = N pk (R, x, v, R' , S) dRdxdv S(x) dx. (B7) 

This result is well understood. As shown in Appendix A. 3, variable v can only take the same value as v' and, hence, 
<7q v can only depend on R' so that x = — (<72 dR(a' v) must be null. Finally, following the same procedure 
we can infer the density of peaks for the remaining infinite series of variables V4, 1)5, ... in infinitesimal ranges. But 
this is actually not necessary since the integration of equation (B4) over these latter variables not entering in P is 
trivial, arriving to 

N^ st (R-^R', x, v, x, S^Sf) dRdR'dxdvdx = N pk (R, x, v, x, R' , 5) P(R', 5 f \v f ) dRdR'dxdPdx, (B8) 

with Vf — v — xAS/(xa' Q ) and P(R',Sf\Pf) given, for variables with subindex /, by equation (A27). Then, by 
integrating equation (B8) over x and v, which leads to 

Np k st {R-> R',x,S -» S f )dRdR' dx = M ^ N(R', S f ) N pk (R, x, S\R' , Sf) dR', (B9) 

taking the Taylor series expansions around S, keeping first order terms in AS, and integrating over x (in the positive 
range) we obtain (see eq. [A13]) 

Np k st (R ->R',S-> Sf) dRdR' = N pk (R, S) P(R', S\R, S) f 1 - d Sf ln[N(R', S f ) N pk (R, S\R' , Sf)] AsldRdR'. 

I Sf—S J 

(BIO) 

Equation (BIO) can also be written, to first order in AS, as 

N^ st (R^ R',6 -» 6 f ) dRdR' = N pk (R,S)dRP(R',S f \R,S)dR', (Bll) 
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with 



P(R',5 f \R,5)dR> = ^N(R\8 f )dR' Np i R ^ R 'j. 6f) (B12) 



giving the probability (only for R' > R; the case R' — R being excluded) that a (non-fissionable) peak with 5 on 
scale R is located on a disjoint background with Sf on scales between R' and R' + dR 1 . This is just what one would 
expect from the same arguments leading to the nesting probability (A13). Notice however that, in contrast to that 
case, in which Sf = S so that the corresponding probability could only include the nesting effect, the probability 
on the left hand side of equation (B12) corresponding to Sf ^ S will now include not only the nesting of peaks 
having evolved from R to Rf < R' but also direct evolution of peaks from R to R' . While, by construction, the 
right-hand member of equation (BIO) only includes the former kind of effect. However, as readily seen from a 
similar development as that leading to N^ st , the density of peaks with R at S which directly evolve into the scale 
R' at Sf is of high order in AS/(R' — R). Therefore, these two expressions coincide, for a given fixed difference 
R' — R, to first order in AS as used in equation (B12). (This reasoning shows, in particular, that eq. [Bll] is not 
true in general; it is only approximately satisfied for very small values of AS. The reason for this is clear. The 
spatial shift of evolving peaks makes the different order scale derivatives of the density contrast of the evolved peak 
at Sf deviate from those of the initial point at the same Sf for arbitrarily large values of AS (see eq. [B3]). The 
equality only holds to first order in AS. Therefore, the probability to find a peak with R at S located in a disjoint 
background with R! at Sf is different, in general, from the probability that the corresponding evolved peak at Sf 
is located in that background. These two probabilities only approximately coincide for very small AS.) 

We are now ready to calculate the net density of peaks with 5 on scales between R and R + dR becoming 
nested into non-nested peaks with Sf = S — AS on scales between R' and R' + dR' , N d (R — > R' , S — > Sf) dRdR '. 
(Superindex d stands for destruction since these filtering events correspond to true mergers, hence the destruction, 
of objects.) To do this we must compute the density of peaks which, after evolving (without fissioning) from S to 
S f , turn out to be nested into non-nested peaks with such larger scales minus the density of (non-fissionable) peaks 
at S which are nested in the ancestors of those disjoint backgounds at Sf. (These ancestors concern, in principle, 
any kind of evolution, although the contribution of direct evolution can be ignored in practice since being of higher 
order in AS). This correction can be readily performed following the same procedure as leading to equation (1). 
The result is the integral equation 

N d (R ->R',5-> S f ) = N^ st (R ^R',S^ S f ) f dR" M ( R "^ N d (R" -» R',S^ S f ) N pk (R, S\R", S) (B13) 

JR P 

for the unknown function N d (R -> R',S-> 5 f ) in terms of N p f st (R -> R',5^> S f ) (eq. [BIO]) and N pk (R, S\R", S) 
(see eq. [1]). 

But we do not need to solve equation (B13). Given the meaning of N d and N, we have N d (R" — > R',5 — ► 
S) = N(R',5)S(R' - R"). (Notice that R" reaches the value R' inside the integral of eq. [B13].) And since 
N pk St ( R -> R',8^> 5) is equal to N pk {R, S) P(R', S\R, S) (see eq. [BIO]) we are led to N d (R -> R', 5 -> S) = for 
R < R' as in the present case. Thus, equation (B13) reduces to a simple relation between first order terms in AS. 
By dividing this relation by AS we are led (see eqs. [BIO] and [A13]) to the Volterra type integral equation of the 
second kind 

N d (R -» R' ,S) = - d Sf [N(R' ,8 f ) N pk (R,S\R' ,S f )\ - /* dR" Mi^D. ^(fl" R' , S) N pk (R, S\R" , S) 

P s f= s JR P 

(B14) 

whose solution gives the wanted net density N d (i? — ► R',S)dRdR'dS of peaks at 6 with scales between 
R and R + dR becoming nested (merging) into non-nested peaks with scales between R' and R' + dR' in 
the next — dS. Equation (B14) can be solved numerically by iteration from the initial approximate solution 
— M(R') p~ x ds f [N p k(R', S) N(R, S\R', 5)]\ g _ g . Actually, this latter function is a very good approximation to the 

wanted solution (with an error of less than 10% for the whole range of R' in the least favorable case of R close to 



Growth History of Objects 



23 



the lower limit of validity) for all power spectra analyzed. Thus, we can simply take 

M{R') 



N d (R->R',S) 



d Sf [N(R',S f )N pk (R,S\R',S f )} 



5 f =6 



(B15) 



with dgN(R', 5) given by numerical solution of the new Volterra type integral equation 

„, M(R") 



d s N(R',S) 



d s N pk (R',5) 



dR" 



■N{R",S)d s N pk (R',S\R",5) 



L 



IR> P 

resulting from differentiation of equation (1). 



dR" M(Jr) d s N(R", 5) N p k(R', 6\R", 6) 



(B16) 



APPENDIX C: NET DENSITY OF NON-NESTED APPEARING PEAKS 
The net density of peaks appearing at (5 with scales between R and R+dR and variables x, V2,v 3l ... in infinitesimal 
ranges is equal to the density of (non-fissionable) peaks with these characteristics minus the density of peaks at 5 
with identical characteristics arising by continuous evolution (hence, without fissioning) from peaks at Si = S + AS, 
with AS positive and arbitrarily small. Therefore, the net density of non-nested peaks appearing (forming) at S 
with scales between R and R + dR is 



N^(R, (5, x, V2, v 3 , •••) AS dRdxdv2 dv 3 .... = 1 — P(R, S, x, v 2 ,v 3 , ... 



x N pk (R, x, V2, v 3 , 5) — N pk (R, x, w 2 , v 3 , 5) dRdxdv2dv 3 



(CI) 



with P(R,S,x,v 2 ,v 3 , ...) the integral (over R 1 > R) of the differential nesting probability given by equation (A14) 
and N^(R, x,V2, V3, S) dRdxdv2 dv% ... the density of evolved peaks, equal to 



Npl{R,x, V2,v 3 ,...,S) dRdxdv 2 dv 3 ... = N pk (R i ,x i ,V2 i ,v 3i , ...,Si) dRi dx t dv 2i dv 3i ... 



(C2) 



(Supcrindex / stands for formation since these filtering events correspond to the formation of new objects.) Equation 
(C2) is but equation (Bl) for the present notation, with the relations between variables with and without subindex 
i equal to the inverse of those given in equations (B2) and (B3). Note that, for the same reasons as in Appendix 
B, we are forced to follow the evolution of the whole infinite set of variables x, v 2 , v 3 , ■■■■ 

From equations (Cf) and (C2) we have 



N^(R, S, x, V2, v 3 , ...) AS dRdxdv2 dv 3 ... = I — P(R, S, x, V2, v 3 , ...) 



N pk (R,x,v 2 ,v 3 , ...,S) - N pk (R i ,x i ,V2i, —,5i)\J\ 



dR dx dvo dv-i . . . 



(C3) 



with J the Jacobian of the transformation from variables with subindex i to variables without subindex. Thus, by 
integrating equation (C3) over all intermediate variables we obtain 

N' (R, S) AS dR = I^N{R, S) - J [1 - P(R, S, x, v 2 , ...)] N pk {R l ,x l ,v 2l , S l ) 



XiRa 2 \ o. 



dxi dv2i 



dR. 



(C4) 



In deriving equation (C4) we have used the fact that N pk (R, S) times the probability that a peak with S on scale 
R be non-nested, 1 — P(R, S), is just equal to N(R, S) (eqs. [Af 3] and [1]). We have also transformed the variables 
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at S in the integral on the right hand side to variables at Si, which balances the Jacobian |J|, and then changed 
the differential dRi by dR according to equation (B2). 

For identical reasons as in Appendix B, it is convenient to express N pk and the nesting probability P in equation 
(C4) in terms of variables v and x instead of v 2 and V3. After this substitution, the integration over variables 
V4,v 5 , ... becomes trivial, and equation (C4) can be written as 

N f (R,S)ASdR = \n(R,S) - ( [1 - P(R, 8, v)] N pk {R h x h v h x h 5 t ) 1 ^— ( R^ - dx l dv l dx\dR, 

L J l XiRa 2 \ a 2 RJ\ J 

(C5) 

where Ri and Si are functions of R and Xi (see cq. [B2]) and <5, respectively. The integral involving only the density 
Npk is straightforward, whereas that involving the product P N p k can be calculated following the same steps as in 
Appendix B. By integrating over the variables Xi, Pi we are led to 



N f (R, 5) ASdR = j A(i?, 8) - J [1 - P(5\5i, x u Ri)] N pk {R t , x u 8 t 



AS ( R ol I 



XiRa 2 V °2 R 



dx t \ dR, (C6) 



with P(Sf\5, x, R) dR the integral over R' of the probability P(R', S\5i, Xi, Ri) dR that a (non-fissionable) peak with 
Ri and Xi at Si is located in a disjoint background with S on scales between R' and R' + dR' (eq. [B12] but with 
the extra dependence on Xi through the factor N p k(R i ,5i,Xi\R,8)/N p i : (Ri,5i,Xi)). Then taking the Taylor series 
expansions around 5, keeping first order terms in A5, integrating over x t (in the positive range), and finally dividing 
by — AS and taking the limit AS — > we arrive to the relation 



<5,=<5 



Nf(R,S)d5dR= \d s N pk (R,S)- f dR'^-^N{R',S)d St N pk {R,8 l \R',S) 

I J R P 

-d R Af pk (S,R)+ f dR' ^^N(R',S)d R Af pk (S,R\S,R')Xd5dR. (C7) 

JR P J 

For simplicity, we will take 

N f (R,S)d5dR^ [d s N pk (R,S) - d R Af pk (S,R)} dSdR (C8) 

which is a reasonable approximation (with an error of less than 15 % and 20 % for scales within the range of validity 
of the model in the cases of the CDM and the n = —2 power law spectra, respectively) to the exact relation (C7). 
Note that, given the meaning of the densities N pk (R, 5) dR and Af pk (S, R) dS (see paper I), expression (C8) is just 
the density of peak trajectories leaving minus entering in the square infinitesimal area of the S vs. R diagram, 
equal to the density of peak trajectories appearing minus disappearing (i.e., fusioning minus fissioning) in that 
inifitesimal area. That is, the approximation leading to (C8) is but the neglect of the correction for nesting. 

APPENDIX D: MASS ACCRETION RATE 
From equation (B2) we have that the mass accreted from S to Sf = 5 — AS, with AS positive and arbitrarily 
small, by any non-nested peak with initial scale between R and R + dR and variable x in an infinitesimal range is 

AM=^^—AS. (Dl) 
dR xa 2 R 

The density of such accreting non-nested peaks is 

N a (R,x,5)dRdx = N(R,x,5)dRdx. (D2) 

In writing equation (D2) we have taken into account that the density of non-nested peaks not accreting because 
merging from 8 to 5— AS is a higher order correction (see the discussion leading to eq. [B14]). Therefore, by dividing 
the density (D2) by N(R, S) dR we obtain the conditional probability p(x, R\5, R) dx that an accreting non-nested 
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peak at S with scale R has the appropriate value of x in order to increase its mass by AM given by equation 
(Dl) in the passage from S to 5 — AS. From equation (1) but for peaks with variable x in an infinitesimal range 
(or, equivalenly, from eq. [A14], after integrating the product N pk P over the remaining variables) this conditional 
probability takes the form 



p(x, R\S, R) dx 



dx 



N(R,S) 



N pk {R,x,S) 



R P 



(D3) 



with N pk (R,x,S) — Af pk (5,x, R) x <r 2 R and N pk (R,x,S\R',S) = J\[ pk {5,x, R\S, R') x <r 2 R in terms of the analogous 
density functions calculated by BBKS (see eq. [A8] and paper I). 

By changing variable x into M' — M + AM we can compute the instantaneous accretion rate of objects of 
mass M per specific range of mass M' of the final object similar to the merger rate (6). By doing so we arrive 
to the fact that this transition rate is identically null. The reason for this is that, as mentioned in Appendix B, 
the density of peaks with scales between R and R + dR at 5 which directly evolve into peaks with scales between 
R' and R' + dR' at Sf is of higher order than one in AS = Sf — S. This result is well understood. From the 
viewpoint of the accreting object, the process is not a transition between two different masses but as a continuous 
mass increase. Consequently, no discrete increment AM can be achieved in the limit AS — > 0. Of course, such a 
continuous evolution of the accreting object during the small interval At necessarily causes (is made at the expense 
of) the merger (capture) of a number of tiny objects which do make a finite transition in mass. And it is taking the 
limit for vanishing At of the change in the number density of these latter objects that one obtains a non-vanishing 
accretion(+merger) rate (see § 3). 

In any event, there is no problem, even from the viewpoint of accreting objects, in obtaining the instantaneous 
mass accretion rate at t for objects of mass M. Equation (Dl) tells us that the instantaneous mass increase rate, 
by accretion, for objects arising from peaks with 5 and the specific value of x is 



dM _dM 1 
dt dR xa 2 R 



dSc 
dt 



(D4) 



Therefore, the instantaneous mass accretion rate at t for objects of mass M (disregarding any other particularity) 
is the average of the specific rate (D4) for the probability function (D3) with R and S = S c expressed in terms of 
M and t through equations (3) and (4). By performing this average we arrive to 



, . dM 
,(M,t) = — 



1 



dR x e f f (72 R 



dS r 



dt 



with 



x eff 



1 



I' OO 

1+ / dR 
Jr 



, M(R') 



N(R',S) U- { S) 



N pk (R,S\R',S) 
N(R,S) 



(D5) 



(D6) 



Since l/(x) is a very good approximation to l/x e ff (with an error of less than 5% for masses within the range of 
validity of the model) we will simply take 



/ , r N dM 



dR (x)a 2 R 



dS c 



dt 



(D7) 
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